Local shell-to-shell energy transfer via nonlocal Interactions in fluid turbulence 
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In this paper we analytically compute the strength of nonlinear interactions in a triad, and the 
energy exchanges between wavenumber shells in incompressible fluid turbulence. The computation 
has been done using first-order perturbative field theory. In three dimension, magnitude of triad 
interactions is large for nonlocal triads, and small for local triads. However, the shell-to-shell 
energy transfer rate is found to be local and forward. This result is due to the fact that the 
nonlocal triads occupy much less Fourier space volume than the local ones. The analytical results 
on three-dimensional shell-to-shell energy transfer match with their numerical counterparts. In 
two-dimensional turbulence, the energy transfer rates to the near-by shells are forward, but to the 
distant shells are backward; the cumulative effect is an inverse cascade of energy. 
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I. INTRODUCTION 

Many equations in physics, e.g., Schrodinger equation and diffusion equation, are local in real space. Here, to 
time-advance a variable at a point, we need the values of the variables and their finite-order derivatives at the same 
point. It is well known that incompressible Navier-Stokes (NS) equation, 

^ + (u-V)u = ~Vp + ^V 2 u, (1) 
V-u = 0, (2) 

is nonlocal in real space 0, Q ■ Here u and p are the velocity and pressure fields respectively, and v is the kinematic 
viscosity. The nonlocality is due to the pressure term of Eq. (TIJ, which is obtained by taking the divergence of 
incompressible Navier-Stokes (NS) equation Q 

V 2 p = -V • {u • Vu} . (3) 

Hence, 

^.o=-/ v '" {u( ty , " > - <4> 

which is nonlocal because p(x, t) depends on the velocity field at x' ^ x. 
In Fourier space, incompressible NS equation is 

dui(k) 2 /i \ * D ... [ dp 

— — Yvk Ui{k) = --F ijm {k) J d u , 



vk Ui(k) = --P ljm (k) J £j3d«j(p)um(q), (5) 
hu t {k) = 0, (6) 



where 



Pijrn (k) — kj ^Sim k 2 ^ km ^ 

k = p + q, and d is the space dimensionality. Note that the factor —kikjk m /k 2 of Py m (k) is due to the pressure term. 
To determine u(k, t+dt) we need the values of field iti(p) where k — p could be quite large. Hence, incompressible NS 
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equation in nonlocal in Fourier space also |2J, |3|. The basic unit of nonlinear interactions in turbulence, called triad 
interactions, involve three vectors (k, p, q) with k = p + q. These triad interactions yield energy transfers among the 
participating modes. The strength of triad interaction is measured using energy exchanges in the triad. 

The energy transfers between two wavenumber shells can be computed using the triad interactions. One of the key 
ingredient of Kolmogorov's phenomenology of turbulence is "local" shell-to-shell energy transfer. That is, maximum 
energy is transferred from a wavenumber shell to the next wavenumber shell. This observation has been verified in 
numerical simulations 0,0]. This result is surprising in view of nonlocal interactions in both real and Fourier space. 
In this paper, we compute the strength of triad interaction, and shell-to-shell energy transfer in incompressible fluid 
turbulence using field-theoretic technique. 

There have been many attempts in the past to compute the strength of triad interactions and the energy transfers 
in fluid turbulence. Kraichnan [f| computed these quantities in both 2D and 3D turbulence using "almost Markovian 
Galilean invariant" turbulence model. He showed the in 3D 35% of the total energy transfer across a unit wavenumber 
sphere involves triads in which the smallest wavenumber is more than one-half of the middle wavenumber. Hence, shell- 
to-shell energy transfer in 3D turbulence is local in wavenumber space. Later Domaradzki and Rogallo numerically 
computed the above quantities and observed that energy transfers in shells are always local, but the triad interactions 
are nonlocal, i.e., triads having three wavenumbers of very different magnitudes have large magnitudes. This is 
succinctly described by Domaradzki and Rogallo as "nonlocal interactions and local energy transfer". Domaradzki 
and Rogallo found their numerical results to be in an excellent agreement with their own Eddy-damped quasi-normal 
Markovian (EDQNM) calculation. They conjectured that the observed energy transfer as being caused by triads 
with at least one wavenumber in the energy-containing range. Ohkitani and Kida Q analyzed the triad interactions 
carefully and concluded that the nonlocal interaction is strong, but the energy exchange occurs predominantly between 
comparable scales. They claimed that the third mode of much larger scale is indifferent to the energy transfer as 
if it were a catalizer in a chemical reaction. Zhou Q numerically computed the energy transfers using different 
wavenumber summation scheme, and found the energy transfers to be local. Waleffe Q did a similar analysis using 
a decomposition of the velocity field in terms of helical modes. Kishida et al. j§| used wavelet basis to address the 
same problem and obtained similar results. For review on this topic, refer to Zhou and Speziale In the present 
paper we re-look at some of the above conjectures. 

In all the above papers, turbulent interactions are measured using a function S(k\p, q) (usually called transfer 
function) that denotes the sum of energy transfers from mode p and q to mode k [IjJ. Dar et al. [12,03 P om t out 
that the energy transfer from one shell to another shell can not be accurately computed using S(k\p,q), essentially 
because the third mode of the interaction could lie outside both the shells under consideration. To overcome this 
difficulty, Dar et al. ^3] modified the above formulation. They used a new function S(k\p\q), called mode-to-mode 
energy transfer rate, for the energy transfer from mode p to mode k, with mode q acting as a mediator, and showed 
that the shell-to-shell energy transfer can be correctly computed by this formalism. The computation of S(kMg) is 
done either numerically or using field-theoretic methods. Kishida et al. used similar formalism as Dar et al. |"12l | for 
wavelets and numerically computed the shell-to-shell energy transfer rates. In the present paper we quantify the triad 
interactions using S(k\p\q), and compute them using field-theoretic method. We also calculate the energy transfer 
rates between wavenumber shells using a first-order perturbation theory. Our analytic arguments justify Domaradzki 
and Rogallo's Q, Zhou's and Ohkitani and Kida's Q numerical results that the turbulent interaction is nonlocal, 
but the shell-to-shell energy transfer is local. 

In this paper we also compute mode-to-mode energy transfer rate S(k\p\q) for space dimension other than 2. For 
d = 2, S(k\p\q) < for most of p < k, unlike 3D case. This property of S(k\p\q) is the reason for the inverse energy 
cascade. We find that the transition from backward to forward energy transfer takes place at d c — 2.25. 

The organization of our paper is as follows: in Sec. 2 we compute S(k\p\q) in the inertial range using first-order 
field theory. It is shown that nonlinear interactions in incompressible NS are nonlocal. The nature of S(k\p\q) for 2D 
and 3D are contrasted. Sec. 3 contains estimates of the shell-to-shell energy transfer for neighbouring and distant 
shells; since the maximal energy transfer takes place between neighbouring shells, the shell-to-shell energy transfer is 
said to be local. In Sees. 4 and 5 we compute shell-to-shell energy transfer rates in 3D and 2D respectively. Section 4 
also contains a comparison of analytical results with their numerical counterparts. Section 6 contains an elementary 
discussion on the energy transfer rates in Burgers turbulence. Section 7 contains conclusions. 

II. NONLOCAL INTERACTIONS IN INCOMPRESSIBLE FLUID TURBULENCE 

Kraichnan 0] has computed magnitudes of triad interactions using transfer function S(k'\p, q) (k' + p + q = 0). 
In the following discussion we will compute the strength of triad interaction using Dar et al.'s mode-to-mode energy 
transfer rate S(k'\p\q) [12] 

5(k'|p|q) = -3 ([k' • u(q)] [u(k') • u(p)]) (7) 
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Figure 1: The interacting triad (k, p, q)/fc = (l,v,w) under the condition k = p + q is represented by a point (v, w) in the 
hatched region. The axis (v',w') are inclined to axis (v,w) by 45 degrees. Note that the local wavenumbers are ti«l,ii)Wl 
or v' w w' « l/y/2. 



that represents energy transfer mode from p to mode k, with mode q acting as a mediator. Here 3 represents 
the imaginary part of the argument. Here we compute the ensemble average of S, (5(k'|p|q)), using the standard 
field-theoretic technique [ill O ITr|. We expand the (5(k'|p|q)) [Eq. 0] to first order in perturbation (see Verma 
0HG3] for details). We assume the flow to be homogeneous and isotropic. 

Following the standard field-theoretic procedure, we perform average of 5(k'|p|q) and obtain an expression for the 
energy transfer rate from mode p to mode k, with mode q as a mediator 

(S(k'\p\q))) = Tl{k ' P ' ^ C ^ C ^ + T2 ( k > p > g) g ( fc ) g (g) + T ^P' <l)C{k)C{ P ) 

v{k)k 2 + v{p)p 2 + v{q)q 2 

where C(k) is the equal-time correlation function, and v(k) is the effective viscosity. The functions Ti(k,p,q) are 
given by 

Ti(k,p,q) = kp((d-3)z + (d-2)xy + 2z 3 + 2xyz 2 + x 2 z) , (9) 
T 2 (k,p,q) = -kp((d-3)z + {d-2)xy + 2z 3 + 2xyz 2 +y 2 z) , (10) 
T 3 (k,p,q) = -kq(xz -2xy 2 z - yz 2 ) , (11) 

where d is the space dimensionality, and x, y, z are the cosines of angles between (p, q), (k, q), and (k, p) respectively. 
We take Kolmogorov's spectrum for the correlation function, i. e., 

C{k} ~ S d (d-1) k^ ' (12) 

and renormalized viscosity of McComb and Watt [3] for v{k) to be 

v{k) = V^^|n| 1/3 fc- 4 / 3 , (13) 

where Kro is Kolmogorov's constant, II is the energy flux, and v* is a constant related to the renormalized viscosity 
0,0]. Note that in 2D fluid turbulence, n is negative for the wavenumber region with 5/3 spectral index. McComb 
and Watt [T5 Hfl| have computed v* using renormalization technique. Here we take Kk = 1.6, v* = 0.38 for 3D, 
and K Ko = 6.3, v* = -0.6 for 2D. 

The interactions are self-similar in the inertial range, which is the region of our interest. Therefore, it is sufficient 
to analyze S(k'\p\q) for triangles (l,p/k, q/k) = (1, v, w). Since, \k — p\ < q < k + p, |1 — vj < w < 1 + v, hence any 
interacting triad (1, v, w) is represented by a point (v, w) in the hatched region of Fig. [ll| . 

The local wavenumbers are v w l,w ~ 1, while the rest are called nonlocal wavenumbers. We substitute C(k) and 
v(k) in Eq. JSJ, which yields 



(S(v,w)) = 



4(27r) 2d ^ 2 H 



S 2 {d- l) 2 k 2d v 

ti(v,w)(vw)- 5 / 3 -( d -V +t 2 {v 1 w)w- 5 / 3 -^-^ +t 3 (v,w)v- 5 / 3 -( d -V 



2/3 _l ,,,2/3 ( 14 ) 
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where ti(v, w) = Ti{k,p, q)/k 2 . For convenience, (S(v', w')} are represented in terms of new variables (v 1 , w') measured 
from the rotated axis shown in the figure ^ It is easy to show that v = 1 + (v' — w')/V2, w = [v 1 + w')/^/2. 

Fig. |21 illustrates density plots of (S(v', w')) without the bracketed factor. Fig. (a) shows the plot for 3D, while Fig. 
(b) shows the one for 2D. Note that (S(v'(v, w), w'(v, w))) is the energy transferred from mode p = v to mode k = l. 
In the white region (positive), energy is transferred from mode p to mode k, while in the dark regions (negative), 
mode p receives energy from mode k. The value of S at (v,w) = (1,1), or (v',w r ) — (1/V2, l/v2) is zero in both 2D 
and 3D. 

In 3D the triads with v 0,w — ► 1 (v' » 0, w' » v2: the top-left corner in v'-w' plot) have large and positive (S), 
implying that the large wavelength modes give a large amount of energy to the modes near k « 1. These observations 
prove that the nonlinear interactions in incompressible NS are nonlocal in Fourier space. In 2D, the triad with v » 
have large-negative (S) implying that the large wavelength modes take energy from the modes near k ~ 1. These 
observations indicate that the interactions in 2D turbulence are nonlocal as well, but the large wavelength modes are 
the sink of energy. 

Another common behaviour in both the dimensions is for w — > 0,v ~ 1 (v' w 0, w' » 0, bottom- left corner in v'-w' 
plot). Here (S) ^> for w < 1, but (S) <C for v > 1, This implies that for these types of triads, the p modes with 
magnitudes less than k always give energy to the k modes, while the p modes with p > k always take energy from the 
k modes. When p, q are much larger than k (v, w — + oo), (S) is small, implying that they interact weakly with k w 1. 
In 3D, (S) for most of these modes are negative implying that they receive energy from k — 1. In 2D, however, (S) 
for a large fraction of these triads are positive; hence they supply energy to k — 1. 

The energy cascade is backward in 2D. This is due to the above mentioned backward energy transfer from k = 1 
mode to the smaller wavenumber modes ((S(v w 0, w w 1)) < 0), and backward energy transfer from large v, w modes 
to k = 1 mode ({S(v,w^$> 1)) > 0). It is interesting to contrast this behaviour with 3D case where (S(v, w)) is 
somewhat opposite to 2D case. 

The function (S(v,w)) in the region with v — * is primarily positive for d = 3, but is negative for d = 2. The 
transition of negative (S) to positive (S) for the region with v — * occurs near d c = 2.25. Please refer to Fig. |^a) 
for the illustration. It can be shown using field-theoretic calculation that the renormalized viscosity vanishes near 
d c = 2.25, and the direction of energy cascade changes from negative to positive at d = d c . Fournier and Frisch (2(1 | 
report d c = 2.05 which differs a bit from our d c - The difference could be because of the fact that Fournier and Frisch 
[2£| use combined energy transfer S(k\p, q) for their EDQNM calculation. In Fig. OJb) we also show the density plot 
of (S(v, w)) for d — 100. This is a representative illustration for large space dimensions. 

The function (S(v, w)) can be estimated in the limiting cases using the method given in Appendix of Leslie In 
the nonlocal region with v — * (naturally, tu»l), 

(S(v,w)) oc vi~ d . (15) 

Clearly, S(v, w) — > oo for both 2D and 3D. This observation is consistent with our earlier observation that interactions 
are nonlocal. However when v « w 1, we find that 

lat \\ / ford = 3 , 

(S(v, w)) oc | a [{v 1)2 _ i {w 1)2 _ i {v _ 1)(w _ 1}] for d = 2 (16) 
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(a) (W 



Figure 3: Density plot of {S(v , to')) of Eq. (0 for (a) d = 2.25 and (b) d = 100. 

This result shows that the interactions within the local triads (fl w 1,wrj 1) are weak. In 3D, among these triads, 
the ones with v < 1 or p < k have (S) > 0, hence the energy transfer is from mode p to mode k; the sign of (S), 
consequently the sign of energy transfer, is reversed for the modes with v > 1. However in 2D, S is somewhat complex 
in the neighbourhood of v = W = 1. 
For v ~ 1 and w — > 0, 

(£(«,«;)) oc (17) 

which again diverges. When v,w — ► oo, 

(5(w,io)) octT^. 

implying that interactions with large wavenumber modes are weak. 

The above estimates are consistent with the graphical plots shown in Figs. and 03 After this discussion, we move 
on to compute the shell-to-shell energy transfer rates in fluid turbulence. 



III. LOCAL SHELL-TO SHELL ENERGY TRANSFER IN INCOMPRESSIBLE FLUID TURBULENCE 

The wavenumber space is divided into shells (koS n , kos n+1 ), where s > 1, and n can take both positive and negative 
values. The energy transfer rate from mth shell (kos m , kos m+1 ) to nth shell (kos n , kos n+1 ) is given by 0] 

T mn = E (S(k\ P \q))- (18) 

If the shell-to-Shell energy transfer rate is maximum for the nearest neighbours, and decreases monotonically with 
the increase of \n— m|, then the shell-to-shell energy transfer is said to be local. 

If the amplitudes of the Fourier modes u(k) are available, either from experiments and or from numerical simulations, 
then we can easily compute the shell-to-shell energy transfer rates using Eqs. f7l lT8|. I n this paper, we compute the 
energy transfer rates between the wavenumbers shells to first order in perturbation |llL IT^ that yields 



T — 



v{k)k 2 + v(p)p 2 + v{q)q 2 
We nondimensionalize the above equations using (Tll | 



dk f dp 

i 

) +TJk n a}n(k)n(n} 

(19) 



k s n <k<k s"+ 1 (27r) d J k s m <p<k s m + 1 (27r) d 
Tifap, q)C{p)C(q) + T 2 {k,p, q)C(k)C{q) + T 3 (k,p, q)C{k)C(p) 



a a a 

k = —: p = — v; q=—w, (20) 
u u u 
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where a is an arbitrary constant wavenumber. For our calculation we choose a = kos n . Three dimensional integral 
under the constraint that k' + p + q = is given by 0] 



/ dp = S d -x ( dpdq (ZL) d 2 (sin af- 3 . (21) 

Jp+q+k=0 J \ K ' 



Using these substitutions, we obtain 



|n| " u {d-iYs d v 

where F(v,w) is given by 



K V \ A 4 ^ 2 'c i t - r ' dv f 1+V dw(vw) d - 2 (sina) d - 3 F(v,w), (22) 



\l-v\ 



tt(v,w)(vw) a d + t 2 (v,w)w => d + t 3 (v,w)v a d 

b{V,w) = j- 575 77-pr: (23) 

V > (l _|_ V 2/Z + y,2/3) V ) 

with ti(v,w) — Ti(k,kv,kw)/k 2 . Eq. II22H provides us with the shell-to-shell energy transfer rates with relative to 
energy flux II. Clearly, T nrn depends only onn-m, or T nm — T n _i iTn _i where i is an integer. Hence, the shell-to-shell 
energy transfer is self-similar. 

Now let us estimate the shell-to-shell energy transfer rates when to <C n. The triads (l,v,w) with v — > participate 
in this energy transfer. As seen in previous section, for these triads S(v,w) cx w^ 4 ^ 3 ^~ d . Therefore, 



d:U 

u 

4/3 



l+v 
\l-v\ 



• d - 2 -- d ' dww d - 2 



where vq is the outer radius of shell m. Since vq — > 0, T nm vanishes. Hence, we are able to show that shell-to-shell 
energy transfer between distant shells in negligible in spite of large interactions between distant wavenumbers. This 
result is essentially due to small volume of v or m-th shell. 

We can also compute the shell-to-shell energy transfer rates for close- by shells. Here the participating triads will 
satisfy v » w « 1. As shown in the earlier section, S(v, w) « — (v — 1) for these triads. Therefore, 

T nm J ~ J d y J dw(vw)(v — 1), (24) 

with n « to. Since our bins are uniform in logarithmic scale, the volume of wavenumber shells is of the order of 1 
when 1. Also the range of v — 1 is of the order of v, which is close to 1. Therefore T nm will be finite. Hence, 

the shell-to-shell energy transfer rates between close-by shells is finite. This result is consistent with the local energy 
transfer assumption of Kolmogorov. 

For close-by shells, p as k, but q can take any value from \k — p\ to k + p. However, it can be easily shown that 
the mode p in the triads with q <C k do not contribute significantly to the shell-to-shell energy transfer. Since, 
S(v,w) cx w'^ for w — * 0, 

Tnm ^ f — f dufi [ dww d - 2 wi~ d w (25) 



u 

1/3 



(26) 



where v = 1 + (iw, and dv = fxdw. Clearly the above integral goes to zero. That is, the shell-to-shell energy transfer 
between close-by shells receive insignificant contribution from the triads satisfying q <C p ~ k. 

In this section we showed using the limiting values of S(k\p\q) that the shell-to-shell energy transfer is local in spite 
of nonlocal interactions among wavenumber modes. However, T nm can be easily computed for the shells in the inertial 
range. In the next section we will do these calculations. 



IV. COMPUTATION OF SHELL-TO-SHELL ENERGY TRANSFER IN 3D 

The shell-to-shell energy transfer rates have been calculated earlier by Ohkitani and Kida 01, and Zhou and Speziale 
[Tnl | using numerical simulation and EDQNM approximation with S(k\p, q). In this section we compute the shell-to- 
shell normalized energy transfer rates T nm /H in three dimensions using Eq. (12:21 . We take s = 2 1 / 4 . The integration 
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Figure 4: Plot of normalized shell-to-shell energy transfer T nm /H vsn-m for d = 3. The nth shell is (k s n , k s n+1 ) with 
s = 2 1//4 . The energy transfer is maximum for n = m± 1, hence the energy transfer is local. The energy transfer is also forward. 



has been done numerically using Gauss-quadrature method. The constants v* — 0.38 and K — 1.6 have been taken 
from McComb and Watt 181 and Verma UJ,UJ- The shell-to-shell energy transfer is self-similar, i. e., T nrn is function 
of n — m. Therefore we compute T nm /H for various n — m. Fig. 0] contains this plot. Note that the shells m and n 
have been assumed to be inside the inertial range. 

From Fig. 0we can infer that the transfer rates T nm in the inertial range are negative for n < m, and positive for 
n > m. Hence a shell gains energy from the smaller wavenumber shells, and loses energy to the higher wavenumber 
shells. This means that the energy cascades from the smaller wavenumbers to the higher wavenumbers (forward). The 
most significant energy transfer takes place from rotom+1. Hence, the shell-to-shell energy transfer is forward and 
local, consistent with Kolmogorov's picture of turbulence. Note that the energy transfer is local in spite of nonlocal 
triad interactions. 

To validate our theoretical calculations, we have also computed the shell-to-shell energy transfer rates using the data 
from Direct Numerical Simulation on a 512 3 grid. The computation was performed when the turbulence was well devel- 
oped. The Reynold's number using Taylor's microscale was 64.8, and skewness was —0.54. We divide the wavenumber 
space into 15 shells with boundaries at wavenumbers (2, 4, 8, 11.3, 13.5, 16, 19, 22.6, 26.9, 32, 38.1, 64, 76.1, 108, 128, 256). 
In the inertial range (k rs 10 — 35) , the shell boundaries are k n — 2(" +11 )/ 4 . Please refer to Dar et al. [l2| for details 
on numerical procedure. Fig. shows the plots of T nm /H vs. n — m for m — 4. .9, which are in the inertial range 
shells. The plots show self-similarity, local, and forward energy transfer for the inertial range shells. The numerical 
and theoretical values are in close agreement. 

For thicker shells s = 2 1 / 2 , the ratio of the smallest to largest wavenumber of the triad is 2\/2. For these shells, 
energy transfer to the nearest neighbouring shell is close to 35%. These numbers are consistent with Kraichnan Q 
and Zhou's 1| results. 

In the next section we will discuss shell-to-shell energy transfer in 2D turbulence. 



SHELL-TO-SHELL ENERGY TRANSFER IN 2D FLUID TURBULENCE 



We compute the shell-to-shell energy transfer in 2D following the same procedure as given above. The wavenumber 
range considered is in the inverse cascade regime (E(k) oc fc- 5 / 3 ). We take v* = -0.6 and K = 6.3 0,03]. 

As shown 

in Fig. U3 the energy transfer rates from the shell m to the three neighbouring shells (m+ 1, m + 2, m + 3) are forward, 
and the transfers are negative for all shells n > m + 3. The above result is very similar to Dar et al. [l^l's numerical 
finding on 2D MHD turbulence (Fig. 10 of Dar et al. [jjj). The negative energy transfer from the distant shells are 
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Figure 5: Plots of normalized shell-to-shell energy transfer T nm /Tl vs. n — m for m from 4 to 9. The plots collapse on each 
other indicating self-similarity. 
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Figure 6: Plot of normalized shell-to-shell energy transfer T n m/|n| vs n-m for d — 2 in the inertial range. The energy 
transfer rate from the shell m to the shells m + 1, m + 2, m + 3 is forward, but m + 4 onward it is negative. The net effect of 
all these transfer is the inverse energy flux II. 



due to negative (S(v, w)} for v — > (top-left of Fig. Gib)), and positive (S(v, w)) for v, w ^> 1. The negative (S(v, w)} 
for v — > indicate that large-wavelength modes receive energy from mode k = 1, and positive (<!?(«, «;)) for u, w 3> 1 
indicate that small- wavelength modes give energy to mode k = 1. A careful inspection of Fig. Efb) indicates that the 
forward energy transfer to shells (m + 1, m + 2, m + 3) is due to a narrow region near v = w = 1, or v' = w' = 1/V%, 
where energy transfer is from lower wavenumber to higher number. 

The above results on shell-to-shell energy transfer is consistent with the energy flux picture. Note that 

oo 

11= ^2 (n-m)T nm . 

n=m-\-l 

When we perform the Y^= m +i( n ~ m )(T n m/\Tl\), we obtain —1, consistent with the inverse cascade of energy in 2D 
turbulence. 
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To summarize, in 2D fluid turbulence, the shell-to-shell energy transfer to the neighbouring shells in forward, but 
the energy transfer is backward for the distant shells. The above behaviour is due to forward local transfer and 
backward nonlocal transfers described in Section 2. 

In the next section we will contrast the energy transfers in incompressible fluid turbulence with relative to Burgers 
turbulence (compressible limit). 

VI. LOCALITY ISSUES IN BURGERS TURBULENCE 

As discussed in the introduction, Navier-Stokes equation is nonlocal in real space due to pressure. In Burgers 
equation 

— + (u ■ V)u = zA7 2 u 

at 

the pressure term is dropped with an implicit assumption that the flow velocity is much greater than sound speed, 
or the sound speed is very small. This is the opposite limit of incompressible NS where the sound speed is infinite. 
Hence a very different behaviour is expected for Burgers equation Q. Clearly, to time-advance the velocity field of 
Burgers equation at a point, we need local field, and its first and second derivative. Hence, Burgers equation is local 
in real-space, and nonlocal in Fourier space. 

The formula for mode-to-mode energy transfer S(k'\p\q) (Eq. [Jj) is not applicable for Burgers equation because it 
is compressible (V-u^O) Therefore, the shell-to-shell energy transfer cannot be computed accurately. Note 

however that energy flux can be computed for Burgers equation. The energy flux is multifractal, and H(k) oc fc~ 1//2 . 
Therefore, £^(fc) oc n 2 / 3 fc- 5 / 3 = A:" 2 jlij. 

VII. CONCLUSIONS 

It is known that the nonlinear interactions in incompressible Navier-Stokes equation is nonlocal in real space due 
to the pressure term. In this paper we investigated locality in Fourier space by computing the strength of triad 
interactions using the formula for the mode-to-mode energy transfer. Our calculation is based on first-order field- 
theoretic technique. We take Kolmogorov's 5/3 powerlaw for the energy spectrum, and the renormalized viscosity for 
the effective viscosity. It has been shown that the magnitudes of interactions for the nonlocal triads k « p ^> q and 
k « q^> p are large, while the interactions are small for the local triads k « p fts q. This result shows that nonlinear 
interactions in incompressible fluid turbulence is nonlocal in Fourier space as well. 

The shell-to-shell energy transfer rates have been investigated by many researchers and ourselves. It is a common 
wisdom that the shell-to-shell energy transfer is local, that is, maximum energy transfer takes place between nearest 
shells. We find that local shell-to-shell energy transfer is compatible with the nonlocal triad interactions because the 
local triads occupy more Fourier space volume as compared to nonlocal (k ~ q ^> p) ones. The local shell-to-shell 
energy transfer via nonlocal triad interactions is consistent, as seen by Domaradzki and Rogallo Zhou t 4!J, Ohkitani 
and Kida @, and Zhou and Speziale [HI m their numerical simulations and EDQNM calculations. We have ourselves 
computed shell-to-shell energy transfer numerically; our theoretical results match with numerical results very well. In 
this paper we show this behaviour analytically. The role of the smallest wavenumber mode in the triad is somewhat 
confusing in earlier papers. We have resolved some of these issues. 

We observe interesting behaviour in two dimensions. The shell-to-shell energy transfer rates to the nearby shells are 
forward, whereas the transfer rates to the far off shells are backward. The net effect is a negative energy flux. This 
theoretical result is consistent with Dar et al.'s numerical finding [l^]. The inverse cascade of energy is consistent 
with the backward nonlocal energy transfer in mode-to-mode picture [S'(fc|p|g)]. We also show that the transition 
from backward energy transfer to forward transfer takes place at d c ~ 2.25. 

The contribution of local triads to the effective or renormalized viscosity has been debated in turbulence literature. 
In Yakhot and Orszag's [22j renormalization theory, the renormalized viscosity gets contribution from highly nonlocal 
wavenumber triads. Kraichnan (2^, 0] first raised the above objection, and proposed some alternatives. The local 
energy transfer and and nonlocal triad interaction results discussed in this paper could be of relevance for this issue; 
this aspect needs further investigation. 

To conclude, an application of field-theoretic techniques to turbulence yields interesting results regarding triad 
interactions and shell-to-shell energy transfers. The method described here has also been applied to magnetohydrody- 
namic turbulence, and the results are presented in Verma et al. [25]. Further investigations of locality in compressible 
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turbulence, and other areas of turbulence will provide us useful clues in furthering our understanding of turbulence. 



[1] U. Frisch, Turbulence (Cambridge University Press, Cambridge, 1995). 

[2] L. D. Landau and E. M. Lifsitz, Fluid Mechanics (Pergamon Press, Oxford, 1987). 

[3] S. A. Orszag, in Fluid Dynamics, Les Houches 1973 Summer School of Theoretical Physics (Gordon and Breach, Berlin, 

1977), p. 273. 
[4] Y. Zhou, Phys. Fluids 5, 1092 (1993). 

[5] J. A. Domaradzki and R. S. Rogallo, Phys. Fluids A 2, 413 (1990). 

[6] R. H. Kraichnan, J. Fluid Mech. 47, 525 (1971). 

[7] K. Ohkitani and S. Kida, Phys. Fluids A 4, 794 (1992). 

[8] F. Waleffe, Phys. Fluids A 4, 350 (1992). 

[9] K. Kishida, K. Araki, S. Kishiba, and K. Suzuki, Phys. Rev. Lett. 83, 5487 (1999). 
[10] Y. Zhou and C. G. Speziale, Appl. Mech. Rev. 51, 267 (1998). 

[11] D. C. Leslie, Development in the Theory of Turbulence (Oxford University Press, Claredon, 1973). 
[12] G. Dar, M. K. Verma, and V. Eswaran, Physica D 157, 207 (2001). 
[13] M. K. Verma, Phys. Rep. 401, 229 (2004). 
[14] R. H. Kraichnan, Phys. Fluids 8, 1385 (1965). 

[15] W. D. McComb, The Physics of Fluid Turbulence (Oxford University Press, Claredon, 1990). 
[16] M. K. Verma, Phys. Rev. E 64, 26305 (2001). 
[17] M. K. Verma, Pramana 61, 577 (2003). 

[18] W. D. McComb and A. G. Watt, Phys. Rev. A 46, 4797 (1992). 

[19] M. K. Verma, Phys. Plasmas 8, 3945 (2001). 

[20] J. D. Fournier and U. Frisch, Phys. Rev. A 17, 747 (1979). 

[21] M. K. Verma, Physica A 277, 359 (2000). 

[22] V. Yakhot and S. A. Orszag, J. Sci. Comput. 1, 3 (1986). 

[23] R. H. Kraichnan, Phys. Fluids 30, 1583 (1987). 

[24] R. H. Kraichnan, Phys. Fluids 30, 2400 (1987). 

[25] M. K. Verma, A. Ayyer, and A. V. Chandra, Phys. Plasmas 12, 82307 (2005). 



